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1. Introduction 



I'm honored to be invited as a lecturer at a Les Houches summer school which 
has a great tradition. I remember reading many of the lectures from past summer 
schools when I was a graduate student and learned a lot from them. I'm also 
looking forward to have a good time in the midst of beautiful mountains, even 
though the weather doesn't seem to be cooperating. I'm not sure if I will ever get 
to see Mont Blanc! 

I was asked to give four general lectures on physics beyond the standard 
model. This is in some sense an ill-defined assignment, because it is a vast subject 
for which we know pratically nothing about. It is vast because there are so many 
possibilities and speculations, and a lot of ink and many many pages of paper had 
been devoted to explore it. On the other hand, we know practically nothing about 
it by definition, because if we did, it should be a part of the standard model of 
particle physics already. I will therefore focus more on the motivation why we 
should consider physics beyond the standard model and discuss a few candidates, 
and there is no way I can present all the examples exhaustively. In addition, after 
reviewing the program, I've realized that there are no dedicated lectures on dark 
matter. Since this is a topic where particle physics and cosmology (I believe) 
are likely to come together in the near future, it is relevant to the theme of the 
school "Particle Physics and Cosmology: the Fabric of Spacetime." Therefore I 
will emphasize this connection in some detail. 

Because I try to be pedagogical in lectures, I wiU probably discuss many points 
which some of you already know very well. Given the wide spectrum of back- 
ground you have, I aim at the common denominator. Hopefully I don't end up 
boring you all! 

1.1. Particle Physics and Cosmology 

At the first sight, it seems crazy to talk about particle physics and cosmology to- 
gether. Cosmology is the study of the universe, where the distance scale involved 
is many Gigaparsecs ~ 10^* cm. Particle physics studies the fundamental con- 
stituent of matter, now reaching the distance scale of ^ 10"^'^ cm. How can they 
have anything in common? 
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The answer is the Big Bang. Discovery of Hubble expansion showed that the 
visible universe was much smaller in the past, and the study of cosmic microwave 
background showed the universe was filled with a hot plasma made of photons, 
electrons, and nuclei in thermal equiUbrium. It was hot. As we contemplate 
earlier and earUer epochs of the universe, it was correspondingly smaller and 
hotter. 

On the other hand, the study of small scales d in particle physics translates to 
large momentum due to the uncertainty principle, p ^ h/d. Since large momen- 
tum requires relativity, it also means high energy E ^ cp ^ hc/d. Physics at 
higher energies is relevant for the study of higher temperatures T ^ E/k, which 
was the state of the earlier universe. 

This way. Big Bang connects microscopic physics to macroscopic physics. 
And we have already seen two important examples of this connection. 

Atomic and molecular spectroscopy is based on quantum physics at the atomic 
distance d ^ 10^^ cm. This spectroscopy is central to astronomy to identify the 
chemical composition of faraway stars and galaxies which we never hope to get 
to directly and measure their redshifts to understand their motion including the 
expansion of the space itself. The cosmic microwave background also originates 
from the atomic-scale physics when the universe was as hot as T ~ 4000 K 
and hence was in the plasma state. This is the physics which we believe we 
understand from the laboratory experiments and knowledge of quantum mechan- 
ics and hence we expect to be able to extract interesting information about the 
universe. Ironically, cosmic microwave background also poses a "wall" because 
the universe was opaque and we cannot "see" with photons the state of the uni- 
verse before this point. We have to rely on other kinds of "messengers" to extract 
information about earher epochs of the universe. 

The next example of the micro-macro connection concerns with nuclear phy- 
sics. The stars are powered by nuclear fusion, obviously a topic in nuclear 
physics. This notion is now well tested by the recent fantastic development in 
the study of solar neutrinos, where the core temperature of the Sun is inferred 
from the helioseismology and solar neutrinos which agree at better than a per- 
cent level. Nuclear physics also determines death of a star. Relatively heavy stars 
even end up with nuclear matter, i.e. neutron stars, where the entire star basically 
becomes a few kilometer-scale nucleus. On the other hand, when the universe 
was as hot as MeV (ten billion degrees Kelvin), it was too hot for protons and 
neutrons to be bound in nuclei. One can go through theoretical calculations on 
how the protons and neutrons became bound in light nuclear species, such as deu- 
terium, ^He, "'He, '^Li, based on the laboratory measurements of nuclear fusion 
cross sections, as well as number of neutrino species from LEP (Large Electron 
Positron collider at CERN). This process is called Big-Bang Nucleosynthesis 
(BBN). There is only one remaining free parameter in this calculation: cosmic 
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baryon density. The resulting predictions can be compared to astronomical de- 
terminations of light element abundances by carefully selecting the sites which 
are believed to be not processed by stellar evolutions. There is (in my humble 
opinion) reasonable agreement between the observation and theoretical predic- 
tions (see, e.g., 1 1 1). This agreement gives us confidence that we understand the 
basic history of the universe since it was as hot as MeV. 




Today 

Fig. 1. Possible messengers from early universe. 

We currently do not have messengers from epochs in early universe above the 
MeV temperature. In other words, our understanding of early universe physics 
is not tested well for T > MeV. Yet many of the topics discussed in this school 
are possible messengers from earlier era: dark matter (10'^ GeV?), baryon asym- 
metry of the universe (10^° GeV?), density perturbations (scalar and tensor com- 
ponents) from the inflationary era (10^^ GeV?). These are the energy scales that 
laboratory measurements have not reached to reveal the full particle spectrum and 
their interactions, hence the realm of physics beyond the standard model. Under- 
standing of such early stages of the universe requires the development in particle 
physics, while the universe as a whole may be regarded as a testing ground of 
hypothesized particle physics at high energies beyond the reach of accelerators. 
This way, cosmology and particle physics help and require each other 

1.2. Next Threshold 

There is a strong anticipation in the community that we are just about to reveal a 
new threshold in physics. Let me tell you why from a historical perspective. 

We (physicists) do not witness crossing a new threshold very often, but each 
time it happened, it resulted in a major change in our understanding of Mother 
Nature. 
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Around year 1900, we crossed the threshold of atomic scale. It is impressive 
to recall how much progress chemists have made without knowing the underlying 
dynamics of atoms and molecules. But the empirical understanding of chemistry 
had clear limitation. For example, van der Waals equation of state showed there 
was the distance scale of about 10~^ cm below which the state-of-art scientific 
knowledge of the time could not be applied, namely the size of atoms. Once 
the technology improved to study precision spectroscopy that allowed people to 
probe physics inside the atoms, a revolution followed. It took about three decades 
for quantum mechanics to be fully developed but it forever changed our under- 
standing of nature. The revolution went on well into the 40's when the marriage 
of quantum mechanics and relativity was completed in Quantum ElectroDynam- 
ics. 

Next important threshold was crossed around 1950 when new hadron reso- 
nances and strange particles were discovered, crossing the threshold of the strong 
interaction scale ~ 10"^^ cm. Discovery of a zoo of "elementary particles" led 
to a great deal of confusion for about three decades. It eventually led to the reve- 
lation of non-perturbative dynamics of quantum field theory, namely confinement 
of quarks, dimensional transmutation, and dynamical synraietry breaking of chi- 
ral symmetry. More importantly, it showed a new layer in nature where quarks 
and gluons take over the previous description of subatomic world with protons 
and neutrons. The experimental verification of this theory. Quantum ChromoDy- 
namics, took well into 90's at numerous accelerators PETRA, PEP, TRISTAN, 
LEP, and HERA. 

One more force that is yet to be fully understood is the weak interaction. Its 
scale was known from the time of Fermi back in 1933 when he wrote the first 
theory of nuclear beta decay. The theory contained one dimensionful constant 
Gf « (300 GeV)"2 « (lO"^*" cm)^. Seven decades later, we are just about 
to reach this energy scale in accelerator experiments, at Tevatron and LHC. We 
do not really know what Nature has in store for us, but at least we've known all 
along that this is another important energy scale in physics. If we are not misled, 
this is the energy scale associated with the cosmic superconductor. Just like the 
Meissner effect lets magnetic field penetrate into a superconductor only over a 
finite distance, the cosmic superconductor lets the weak force carried by W and 
Z bosons go over a tiny distance: a billionth of a nanometer. Right now we are 
only speculating what revolution may take place at this distance scale. A new 
layer of matter? New dimensions of space? Quantum dimensions? Maybe string 
theory? We just don't know yet. 

Of course historical perspective does not guarantee that history repeats itself 
in an equally exciting fashion. But from all what we know, there is a good reason 
to think that indeed a new threshold is waiting to be discovered at the TeV energy 
scale, as I will discuss in the next section. Another simple fact is that crossing 
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a new threshold is something like twice-in-a-century experience. I'm excited to 
think that we are just about to witness one, a historic moment. 

An interesting question is what fundamental physics determines these thresh- 
olds. The atomic scale, that looked like a fundamental limitation in understanding 
back in the 19th century, did not turn out to be a fundamental scale at all. It is a 
derived scale from the mass of the electron and the fundamental constants, 

as = — w 10"^ cm. (1.1) 
e^me 

The strong-interaction scale is also a derived energy scale from the coupling con- 
stant 

o, = Me-s^'/'^'^^)^" « 10-13 cm, (1.2) 

where (/, is the strong coupling constant defined at a high-energy scale M and 
bo is the beta function coefficient. Because of the asymptotic freedom, the strong 
coupling constant is weak (what an oxymoron!) at high energies, while it be- 
comes infinitely strong at low energies. The scale of strong interaction is where 
the strength of the interaction blows up. In other words, the two thresholds we 
have crossed so far were extremely exciting, yet they turned out to be not funda- 
mental! They point to yet deeper physics that determine these parameters in na- 
ture. Maybe the weak-interaction scale is also a derived scale from some deeper 
physics at yet shorter distances. 



2. Why Beyond the Standard Model 

2.1. Empirical Reasons 

Until about ten years ago, particle physicists lamented that the standard model 
described every new data that came out from experiments and we didn't have 
a clue what may lie beyond the standard model. Much of the discussions on 
physics beyond the standard model therefore were not based on data, but rather 
on theoretical arguments, primarily philosophical and aesthetic displeasure with 
the standard model. It all changed the last ten years when empirical evidence 
appeared that demonstrated that the standard model is incomplete: 

• Non-baryonic dark matter, 

• Dark energy, 

• Neutrino mass, 

• Nearly scale-invariant, Gaussian, and apparently acausal density perturba- 
tions, 

• Baryon asymmetry. 
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I will discuss strong evidence for non-baryonic dark matter and dark matter later 
in my lectures. Density fluctuation is covered in many other lectures in this 
school by Lev Kofman, Sabino Matarrese, Yannick Mellier, Simon Prunet, and 
Romain Teyssier. Neutrino mass is discussed by Sergio Pastor, and baryon asym- 
metry by Jim Cline. The bottom line is simple: we already know that there must 
be physics beyond the standard model. However, we don't necessarily know the 
energy (or distance) scale for this new physics, nor what form it takes. One con- 
servative approach is to try to accommodate all of these established empirical 
facts into the standard model with minimum particle content: The New Mini- 
mal Standard Model IS). I will discuss some aspects of the model later. But 
theoretical arguments suggest the true model be much bigger, richer, and more 
interesting. 



2.2. Philosophical and Aesthetic Reasons 

What are the theoretical arguments that demand physics beyond the standard 
model? As 1 mentioned already, they are based on somewhat philosophical ar- 
guments and aesthetic desires and not exactly on firm footing. Nonetheless they 
are useful and suggestive, especially because nature did solve some of the similar 
problems in the past by invoking interesting mechanisms. A partial list relevant 
to my lectures here is 

• Hierarchy problem: why Gp 10"^ GeV-^ < Gjv lO'^^ GeV-^? 

• Why Oqcd < < 1? 

• Why are there three generations of particles? 

• Why are the quantum numbers of particles so strange, yet do anomalies cancel 
so non-trivially? 

For an expanded list of the "big questions", see e.g., |[3|. 

To understand what these questions are about, it is useful to remind ourselves 
how the standard model works. It is a gauge theory based on the SU{3) x 
SU{2) X U (1) gauge group with the Lagrangian 

+QiipQi + LiipLi + UiipUi + diipdi + eiipa 
+ (Y^'Q,UjH + Y'^Qd.H + Yi'Ue.H + h.c.) 

+ {D^H)^{Df'H) - \{H^Hf - m^H^H + -J—e^'^P^TiiG^^Gp^). 

oZtt 

(2.1) 
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It looks compact enough that it should fit on a T-shirtOl Why don't we see such a 
T-shirt while we see Maxwell equations a lot? 

The first two lines describe the gauge interactions. The covariant derivatives 
p = ^^D^i_ in the second line are determined by the gauge quantum numbers 
given in this table: 





SU{i) 


5(7(2) 


U{1) 


chirality 


Q 


3 


2 


+1/6 


left 


u 


3 


1 


+2/3 


right 


D 


3 


1 


-1/3 


right 


L 


1 


2 


-1/2 


left 


E 


1 


1 


-1 


right 



This part of the Lagrangian is well tested, especially by the LEP/SLC data in the 
90's. However, the quantum number assignments (especially U{\) hypercharges) 
appear very strange and actually hard to rememberjj Why this peculiar assign- 
ment is one of the things people don't like about the standard model. In addition, 
they are subject to non-trivial anomaly cancellation conditions for SU{2>YU{\), 
SU{2fU{l), gTavity2c/(l), U{lf, and Witten's SU{2) anomalies. Many of us 
are left with the feeling that there must be a deep reason for this baroque quantum 
number assignments which had led to the idea of grand unification. 

The third line of the Lagrangian comes with the generation index i, j — 1,2,3 
and is responsible for masses and mixings of quarks and masses of charged lep- 
tons. The quark part has been tested precisely in this decade at i?-factories while 
there is a glaring omission of neutrino masses and mixings that became estab- 
lished since 1998. In addition, it appears unnecessary for nature to repeat el- 
ementary particles three times. The repetition of generations and the origin of 
mass and mixing patterns remains an unexplained mystery in the standard model. 

The last line is completely untested. The first two terms describe the Higgs 
field and its interaction to the gauge fields and itself. Having not seen the Higgs 
boson so far, it is far from established. The mere presence of the Higgs field poses 
an aesthetic problem. It is the only spinless field in the model, but it is introduced 
for the purpose of doing the most important job in the model. In addition, we 
have not seen any elementary spinless particle in nature! Moreover, the potential 

' It reminds me of an anecdote from when the standai'd model was just about getting off the ground 
around 1978. There was a convergence of the data to the standard model and people got very excited 
about it. Then Tini Veltman gave a talk asking "do you really think this is great model?" and wrote 
down every single term in the Lagrangian without using a compact notation used here over pages and 
pages of transparencies. Unfortunately I don't remember who told me this story. 

-I often told my friends that I chose physics over chemistry or biology because I didn't want to 
memorize anything, but this kind of table casts serious doubt on my choice! 
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needs to be chosen with < to cause the cosmic superconductivity which 
does not give any reason why our universe is in this state. I will discuss more 
problems about it in a few minutes. Overall, this part of the model looks very 
artificial. 

The last term is the so-called 0-term in QCD and violates T and CP. The 
vacuum angle 9 is periodic under 9-^6 + 2tt, and hence a "natural" value of 
9 is believed to be order unity. On the other hand, the most recent experimental 
upper limit on the neutron electric dipole moment |o?„| < 2.9 x lO^^^e cm (90% 
CL) |40 ti-anslates to a stringent upper limit 9 < (1.2 ± 0.6) x 10^^" using the 
formula in [5|. Why 9 is so much smaller than the "natural" value is the strong 
CP problem, and again the standard model does not offer any explanations. 

Now we have more to say about the Higgs sector (the third line). Clearly it is 
very important because ( 1 ) this is the only part of the Standard Model which has a 
dimensionful parameter and hence sets the overall energy scale for the model, and 
(2) it has the effect of causing cosmic superconductivity without explaining its 
microscopic mechanism. For the usual superconductors studied in the laboratory, 
we can use the same Lagrangian, but it is derived from the more fundamental 
theory by Bardeen, Cooper, and Schrieffer. The weak attractive force between 
electrons by the phonon exchange causes electrons to get bound and condense. 
The "Higgs" field is the Cooper pair of electrons. And one can show why it has 
this particular potential. In the standard model, we do not know if Higgs field is 
elementary or if it is made of something else, nor what mechanism causes it to 
have this potential. 

All the puzzles raised here (and more) cry out for a more fundamental theory 
underlying the Standard Model. What history suggests is that the fundamental 
theory lies always at shorter distances than the distance scale of the problem. 
For instance, the equation of state of the ideal gas was found to be a simple 
consequence of the statistical mechanics of free molecules. The van der Waals 
equation, which describes the deviation from the ideal one, was the consequence 
of the finite size of molecules and their interactions. Mendeleev's periodic table 
of chemical elements was understood in terms of the bound electronic states, 
Pauli exclusion principle and spin. The existence of varieties of nuclide was due 
to the composite nature of nuclei made of protons and neutrons. The list could 
go on and on. Indeed, seeking answers at more and more fundamental level is 
the heart of the physical science, namely the reductionist approach. 

The distance scale of the Standard Model is given by the size of the Higgs 
boson condensate v = 250 GeV. In natural units, it gives the distance scale of 
d = hc/v = 0.8 X 10~^^ cm. We therefore would like to study physics at distance 



^This is an amazing limit. If you blow up tlie neutron to tlie size of the Earth, this limit cocesponds 
to a possible displacement of an electron by less than ten microns. 
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scales shorter than this eventually, and try to answer puzzles whose partial list 
was given in the previous section. 

Then the idea must be that we imagine the Standard Model to be valid down 
to a distance scale shorter than d, and then new physics will appear which will 
take over the Standard Model. But applying the Standard Model to a distance 
scale shorter than d poses a serious theoretical problem. In order to make this 
point clear, we first describe a related problem in the classical electromagnetism, 
and then discuss the case of the Standard Model later along the same line 19]. 



2.3. Positron Analogue 

In the classical electromagnetism, the only dynamical degrees of freedom are 
electrons, electric fields, and magnetic fields. When an electron is present in the 
vacuum, there is a Coulomb electric field around it, which has the energy of 

1 

Ai?couiomb = -; ■ (2.2) 

Aneo Te 

Here, Ve is the "size" of the electron introduced to cutoff the divergent Coulomb 
self-energy. Since this Coulomb self-energy is there for every electron, it has to 
be considered to be a part of the electron rest energy. Therefore, the mass of the 
electron receives an additional contribution due to the Coulomb self-energy: 

{meC^)ohs = {'rn^C^)kare + AScoulomb- (2.3) 

Experimentally, we know that the "size" of the electron is small, < 10^^^ cm. 
This impUes that the self-energy /S.E is greater than 10 GeV or so, and hence 
the "bare" electron mass must be negative to obtain the observed mass of the 



electron, with a fine cancellation likq 

0.000511 = (-3.141082 + 3.141593) GeV. (2.4) 

Even setting a conceptual problem with a negative mass electron aside, such 
a fine-cancellation between the "bare" mass of the electron and the Coulomb 
self-energy appears ridiculous. In order for such a cancellation to be absent, we 
conclude that the classical electromagnetism cannot be applied to distance scales 
shorter than / {Aneome,c?) = 2.8 x 10~^^ cm. This is a long distance in the 
present-day particle physics' standard. 

The resolution to this problem came from the discovery of the anti-particle of 
the electron, the positron, or in other words by doubling the degrees of freedom 
in the theory. The Coulomb self-energy discussed above can be depicted by a dia- 
gram Fig. |2] where the electron emits the Coulomb field (a virtual photon) which 



*Do you recognize tt? 
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Fig. 2. The Coulomb self-energy of the electron. 



e; 




Fig. 3. The bubble diagram which shows the fluctuation of the vacuum. 



is absorbed later by the electron (the electron "feels" its own Coulomb field)! 
But now that we know that the positron exists (thanks to Anderson back in 1932), 
and we also know that the world is quantum mechanical, one should think about 
the fluctuation of the "vacuum" where the vacuum produces a pair of an electron 
and a positron out of nothing together with a photon, within the time allowed 
by the energy-time uncertainty principle At ~ h/AE ^ h/{2mf.c^) (Fig. O. 
This is a new phenomenon which didn't exist in the classical electrodynamics, 
and modifies physics below the distance scale d ^ cAt ~ hc/{2mec'^) = 200 x 
10^^"^ cm. Therefore, the classical electrodynamics actually did have a finite ap- 
plicability only down to this distance scale, much earlier than 2.8 x lO^^'^ cm as 
exhibited by the problem of the fine cancellation above. Given this vacuum fluc- 
tuation process, one should also consider a process where the electron sitting in 
the vacuum by chance annihilates with the positron and the photon in the vacuum 
fluctuation, and the electron which used to be a part of the fluctuation remains 
instead as a real electron (Fig.|4]l. V. Weisskopf 1 10| calculated this contribution 
to the electron self-energy, and found that it is negative and cancels the leading 



'The diagrams Figs.|2]|4]are not Feynman diagrams, but diagrams in the old-fashioned perturba- 
tion theory with different T-orderings shown as separate diagrams. The Feynman diagram for the 
self-energy is the same as Fig.|2] but represents the sum of Figs.|2]|4]and hence the linear divergence 
is already cancelled within it. That is why we nomially do not hear/read about linearly divergent 
self-energy diagrams in the context of field theory. 
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Fig. 4. Another contribution to the electron self-energy due to the fluctuation of the vacuum. 



piece in the Coulomb self-energy exactly Q 



1 e 

A-Bpair — — -; 

47reo Te 



(2.5) 



After the linearly divergent piece 1 /r^ is canceled, the leading contribution in the 
^ limit is given by 



Sen f}j 

AE = AiJcoulomb + AiJpair = — mgC^ log . 

47r nif.cre 



(2.6) 



There are two important things to be said about this formula. First, the correction 
/S.E is proportional to the electron mass and hence the total mass is proportional 
to the "bare" mass of the electron. 



1 



3a , h 
— log 

47r rrieCre 



(2.7) 



Therefore, we are talking about the "percentage" of the correction, rather than a 
huge additive constant. Second, the correction depends only logarithmically on 
the "size" of the electron. As a result, the correction is only a 9% increase in 
the mass even for an electron as small as the Planck distance Ve = 1/Mpi — 
1.6 X 10^33 

The fact that the correction is proportional to the "bare" mass is a consequence 
of a new symmetry present in the theory with the antiparticle (the positron): the 
chiral symmetry. In the limit of the exact chiral symmetry, the electron is mass- 
less and the symmetry protects the electron from acquiring a mass from self- 
energy corrections. The finite mass of the electron breaks the chiral symmetry 
explicitly, and because the self-energy correction should vanish in the chiral sym- 
metric limit (zero mass electron), the correction is proportional to the electron 
mass. Therefore, the doubling of the degrees of freedom and the cancellation 
of the power divergences lead to a sensible theory of electron applicable to very 
short distance scales. 



An earlier paper by Weisskopf actually found two contributions to add up. After Furry pointed 
out a sign mistake, he published an enata with no linear divergence. I thank Howie Haber for letting 
me know. 
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2.4. Hierarchy Problem 

In the Standard Model, the Higgs potential is given by 

y = m^li/p + A|i7|^, (2.8) 

where v'^ — {H)"^ = —m?/2\ = (176 GeV)^. Because perturbative unitarity 
requires that A < 1, —m? is of the order of (100 GeV)^. However, the mass 
squared parameter m? of the Higgs doublet receives a quadratically divergent 
contribution from its self-energy corrections. For instance, the process where the 
Higgs doublets splits into a pair of top quarks and come back to the Higgs boson 
gives the self-energy correction 

A"^?op = -6-^^, (2.9) 

where th is the "size" of the Higgs boson, and /ij « 1 is the top quark Yukawa 
coupling. Based on the same argument in the previous section, this makes the 
Standard Model not applicable below the distance scale of 10~^^ cm. This is the 
hierarchy problem. In other words, if we don't solve this problem, we can't even 
talk about physics at much shorter distances without an excessive fine-tuning in 
parameters. 

It is worth pondering if the mother nature may fine-tune. Now that the cos- 
mological constant appears to be fine-tuned at the level of 10~^^°, should we be 
really worried about the fine-tuning of u^/Af|>; « 10~^° f6l? In fact, some peo- 
ple argued that the hierarchy exists because intelligent life cannot exist otherwise 
Q. On the other hand, a different way of varying the hierarchy does seem to 
support stellar burning and life |8|. I don't get into this debate here, but I'd like 
to just point out that a different fine-tuning problem in cosmology, horizon and 
flatness problems, pointed to the theory of inflation, which in turn appears to be 
empirically supported by data. I just hope that proper solutions will be found 
to both of these fine-tuning problems and we will see their manifestations at the 
relevant energy scale, namely TeV. You have to be an optimist to work on big 
problems. 



3. Examples of Physics Beyond the Standard Model 

Given various problems in the standard model discussed in the previous section, 
especially the hierarchy problem, many possible directions of physics beyond the 
standard model have been proposed. I can review only a few of them here given 
the spacetime constraint. But I especially emphasize the aspect of the models 
that leads to a (nearly) stable neutral particle as a good dark matter candidate. 



Physics Beyond the Standard Model and Dark Matter 



17 



3.1. Supersymmetry 

The motivation for supersymmetry is to make the Standard Model applicable 
to much shorter distances so that we can hope that the answers to many of the 
puzzles in the Standard Model can be given by physics at shorter distance scales 
ifTTl . In order to do so, supersymmetry repeats what history did with the positron: 
doubling the degrees of freedom with an explicitly broken new symmetry. Then 
the top quark would have a superpartner, the stop0 whose loop diagram gives 
another contribution to the Higgs boson self energy 

^<op = +6:^^. (3.1) 

The leading pieces in l/r/f cancel between the top and stop contributions, and 
one obtains the correction to be 

H t 

One important difference from the positron case, however, is that the mass 
of the stop, mf, is unknown. In order for the /S.m? to be of the same order of 
magnitude as the tree-level value — —2Xv^, we need to? to be not too far 
above the electroweak scale. TeV stop mass is already a fine tuning at the level of 
a percent. Similar arguments apply to masses of other superpartners that couple 
directly to the Higgs doublet. This is the so-called naturalness constraint on the 
superparticle masses (for more quantitative discussions, see papers [12|). 

Supersymmetry doubles the number of degrees of freedom in the standard 
model. For each fermion (quarks and leptons), you introduce a complex scalar 
field (squarks and sleptons). For each gauge boson, you introduce gaugino, a 
partner Majorana fermion (a fermion field whose anti-particle is itself). I do not 
go into technical aspect of how to write a supersymmetric quantum field theory; 
you should consult some review articles fl3] [T4| . 

One important point related to dark matter is the proton longevity. We know 
from experiments such as SuperKamiokande that proton is very long lived (if not 
immortal). The life time for the decay mode p e+7r" is longer than 1.6 x 
Iq33 years, at least twenty-three orders of magnitude longer than the age of the 
universe! On the other hand, if you write the most general renormalizable theory 
with standard model particle content consistent with supersymmetry, it allows 
for vertices such as eijkU^dPs'^ and eu*s* (here i, j, k are color indices). Then 



'This is a terrible name, which was originally meant to be "scalar top" or "supersymmetric top." 
Some other names are even worse: sup, sstrange, etc. If supersymmetry will be discovered at LHC, 
we should seriously look for better names for the superparticles, maybe after the names of rich donors. 
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one can draw a Feynman diagram like one in Fig. |5] If the couplings are 0(1) 
and superparticles around TeV, one finds the proton lifetime as short as Tp ^ 

mj/rrip ^ 10^^^ sec; a little too short! 



> - < 



Fig. 5. A possible Feynman diagram with supersymmetric particles that can lead to a too-rapid proton 
decay p ^ e+ tt" . 



Because of this embarrassment, we normally introduce a Z2 symmetry called 
"i?-parity" defined by 

Rp = [-lfB+L+2s ^ ^_-^ynattor^^^ (3 3) 

where s is the spin. What it does is to flip the sign of all matter fields (quarks 
and leptons) and perform 27r rotation of space at the same time. In effect, it 
assigns even parity to all particles in the standard model, and odd parity to their 
superpartners. Here is a quick check. For the quarks, B = 1/3, L = 0, and 
s — 1/2, and we find Rp ~ +1, while for squarks the difference lies in s = 
and hence Rp = —1. This symmetry forbids both of the bad vertices in Fig.|5] 

Once the i?-parity is imposedjl there are no baryon- and lepton-number vi- 
olating interaction you can write down in a renormalizable Lagrangian with the 
standard model particle content. This way, the i?-parity makes sure that pro- 
ton is long lived. Then the lightest supersymmetric particle (LSP), with odd 
i?-parity, cannot decay because there are no other states with the same i?-parity 
with smaller mass it can decay into by definition. In most models it also turns out 
to be electrically neutral. Then one can talk about the possibility that the LSP is 
the dark matter of the universe. 



3.2. Composite Higgs 

Another way the hierarchy problem may be solved is by making the Higgs boson 
to actually have a finite size. Then the correction in Eq. (12.91 ) does not require 
tremendous fine-tuning as long as the physical size of the Higgs boson is about 

**An obvious objection is that imposing _R-parity appears ad hoc. Fortunately there are several 
ways for it to emerge from a more fundamental theory. Because the /J-parity is anomaly-free 1151 . 
it may come out from string theory. Or Rp can arise as a subgroup of the 50(10) grand unified 
gauge group because the matter belongs to the spinor representation and Higgs to vector, and hence 
2tt rotation in the gauge group leads precisely to (_i)mattcr jj ^jg^ ^ accidental symmetry 
due to other symmetries of the theory 1161 1171 so that it is slightly broken and dark matter may 
eventually decay. 
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th ~ (TeV)~^ ~ 10^^^ cm. This is possible if the Higgs boson is a composite 
object made of some elementary constituents. 

The original idea along this line is called technicolor (see reviews ITSl [191 ). 
where a new strong gauge force binds fermions and anti-fermions much like 
mesons in the real QCD. Again just like in QCD, fermion anti-fermion pair 
have a condensate {-ipip) ^ breaking chiral symmetry. In technicolor theo- 
ries, this chiral symmetry breaking is nothing but the breaking of the electroweak 
SU{2) X U{1) symmetry to the U (1) QED subgroup. Because the Higgs boson 
is heavy and strongly interacting, it is expected to be too wide to be seen as a 
particle state. 

It is fair to say, however, that the technicolor models suffer from various prob- 
lems. First of all, it is difficult to find a way of generating sufficient masses for 
quarks and leptons, especially the top quark, because you have to rely on higher 
dimension operators of type qqipi/j / A"^ . The scale A must be low enough to gen- 
erate mt, while high enough to avoid excessive flavor-changing neutral current. 
In addition, there is tension with precision electroweak observables. These ob- 
servables are precise enough that they constrain heavy particles coupled to Z- 
and M^-bosons even though we cannot produce them directly]^ 

Because of this issue, there are various other incarnations of composite Higgs 
idea, which try to get a relatively light Higgs boson as a bound state Il26ll27l . One 
of the realistic models is called "little Higgs" [24 25 1. Because of the difficulty 
of achieving Higgs compositeness at the TeV scale, we are better off putting 
off the compositeness scale to about 10 TeV to avoid various phenomenological 
constraints. Then you must wonder if the problem with Eq. (12.9b comes back. 
But there is a way of protecting the scale of Higgs mass much lower than the 
compositeness scale by using symmetries similar to the reason why a pion is 
so much lighter than a proton. If you are clever, you can arrange the structure 
of symmetry such that it eliminates the one-loop correction in Eq. ( 12. 9t and the 
correction arises only at the two-loop level. Then the compositeness ~ 10 TeV 
is not a problem. 

Another attractive idea is to use extra dimensions to generate the Higgs field 
from a gauge field, called "Higgs-gauge unification" ll28l |29l [30l ISTl [32l . We 
know the mass of the gauge boson is forbidden by the gauge invariance. If the 
Higgs field is actually a gauge boson (spin one), but if it is spinning in extra di- 
mensions, we (as observers stuck in four dimensions) perceive it not to spin. Not 
only this gives us raison d'etre of (apparently) spinless degrees of freedom, it 
also provides protection for the Higgs mass and hence solves the hierarchy prob- 
lem. The best implementation of this line of thinking is probably the holographic 

'it is curious that higher dimensional versions of technicolor models called Higgsless models 
1211 do much better 1221 . A supersymmetric version of technicolor also does better than the original 
technicolor |23| . 
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Higgs model in Refs. Il33l [34l which involves the warped extra dimension I will 
briefly discuss in the next section. It should also be said that many of the ideas 
mentioned here are closely related to each other lf35]| . 

Similarly to the case of supersymmetry, people often introduce a Z2 symmetry 
to avoid certain phenomenological embarrassments. In little Higgs theories, tree- 
level exchange of new particles tend to cause tension with precision electroweak 
constraints. Then the new states must be sufficiently heavy so that the hierarchy 
problem is reintroduced. By imposing "T-parity," new particles can only appear 
in loops for low-energy processes and the constraints can be easily avoided [36J. 
Then the lightest T-odd particle (LTP) becomes a candidate for dark matter. In 
technicolor models, the lightest technibaryon is stable (just like proton in QCD) 
and a dark matter candidate [37J. 



3.3. Extra Dimensions 

The source of the hierarchy problem is our thinking that there is physics at much 
shorter distances than 10^^^ cm. What if there isn't? What if physics ends 
at 10~^^ cm where quantum field theory stops applicable and is taken over by 
something more radical such as string theory? Normally we associate the ulti- 
mate limit of field theory with the Planck scale dpi ~ 10"'^^ cm = G]^^ where 
the gravity becomes as strong as other forces and its quantum effects can no 
longer be ignored. How then can the quantum gravity effects enter at a much 
larger distance scale such as 10~^^ cm? 

One way is to contemplate "large" extra dimensions of size R |38|. Imagine 
there are n extra dimensions in addition to our three-dimensional space. If you 
place two test masses at a distance r much shorter than R, the field lines of 
gravity spread out into 3 + n dimensions and the force decreases as the surface 
area . However, if the distance is longer than R, there is a limit to which 

how much the field lines can spread because they are squeezed within the size 
R. Therefore, the force decreases only as r^^ for r ^ R, reproducing the usual 
inverse square law of gravity. It turns out that the inverse square law is tested only 
down to 44 /im |39| (even though this is very impressive!) and extra dimensions 
smaller than that are allowed experimentally. 

Matching two expressions for the gravitational force at r = i?, we can related 
the Newton's constant in ?i + 3 dimensions to the Newton's constant Gat in three 
dimensions 

and hence Gn+3 = GnR^- In the natural unit H ^ c ~ 1, the mass scale 
of gravity is related to the Planck scale by M"+^i?" — Mp^. Even if the true 
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energy scale of quantum gravity is at Af ^ 1 TeV, we may find an apparent scale 
of gravity to be Mpi ^ 10^^ GeV. Then the required size of extra dimensions is 
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(3.5) 



Obviously the n — 1 case is excluded because R is even bigger than lAU ~ 
10^^ cm. The case n = 2 is just excluded by the small-scale gravity experiment, 
while n < 3 is completely allowed. 




Fig. 6. Large extra dimensions. Even though the three-brane is drawn at the ends of extra dimensions, 
it does not have to be. 



If we don't see the extra dimensions directly, what do they do to us? Let us 
look at the case of just one extra dimension with periodic boundary condition 
y ^ y + 2ttR. Then all particles have wave functions on the coordinate y that 
satisfies -0(2/ + 27ri?) = ip{y)- They can of course depend on the usual four- 
dimensional space time x, too. One can expand it in Fourier modes 

V'(x,2/) = ^,An(a;)e^"^/^. (3.6) 

n 

The momentum along the y dkection is py = —idy = n/R, and the total energy 
of the particle is 




(3.7) 



Namely that you find a tower of particles of mass m„ = n/R, called Kaluza- 
Klein states. 
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Of course, the standard model is tested down to 10^ cm, and we have not 
found Kaluza-Klein excitation of electron, etc. This is not a problem if we are 
stuck on a three-dimensional sheet (thiee-brane) embedded inside the n + 3 di- 
mensional space. The branes are important objects in string theory and it is easy 
to get particles with gauge interactions stuck on them. The brane may be freely 
floating inside extra dimensions or may be glued at singularities {e.g., orbifold 
fixed points). The simplest way to use large extra dimensions is to assume that 
only gravity is spread out in extra dimensions, while the standard model particles 
are all on a three-brane. 

Cosmology with large extra dimension is an iffy subject, however The Kaluza- 
Klein excitation of gravitons can be produced in early universe and the cosmol- 
ogy would be different from the standard Friedmann univese (see, e.g., [4Q\ ). 
I will not get into this discussion here. 

Instead of models with large extra dimensions, models with small extra di- 
mensions of size R « 10^^^ cm « TeV^^ are also interesting which allow 
for normal cosmology below TeV temperatures. This would also allow us stan- 
dard model particles to live in extra dimensions, too, because our Kaluza-Klein 
excitations have been too heavy to be produced at accelerators so far There are 
many versions of small extra dimensions. 

One very popular version is warped extra dimension f4T]. Instead of flat met- 
ric in the extra dimensions, it sets up an exponential behavior. It is something 
like Planck scale varies from 10^^ GeV to TeV as you go across the 5th dimen- 
sion. Therefore, physics does end at TeV if you on one end of the 5th dimension, 
while it keeps going to 10^^ GeV on the other end. The hierarchy problem may 
be solved if Higgs resides on (or close to) the "TeV brane." This set up attracted 
a lot of attention because the bulk is actually a slice of anti-de Sitter space which 
has nice features of preserving supersymmetry, leading to AdS/CFT correspon- 
dence 1421 . etc. It is also possible to obtain quite naturally from string theory B3l . 
In a grand unified model from warped extra dimension, the proton longevity is 
an issue which is solved by a Z3 symmetry, and the lightest Z3 -charged particle 
(LZP) is a candidate for dark matter ll44l . 

It is also possible to have the "flat" extra dimension at the TeV scale and put 
all the standard model particles in the 5D bulk, called Universal Extra Dimension 
(UED) ll45l . It is tricky to get chiral fermions in four dimensions if they are 
embedded in higher dimensional space. If you start out with five-dimensional 
Dirac equation 

iij^d^ + j^dy)^ix,y)^0, (3.8) 

'"Historically, unified theories and string theory assumed dpi 10~^^ cm. TeV-sized extra 
dimensions are much larger than this, but I'm calling them "small" for the sake of distinction from 
the large extra dimensions. 
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y=0 y=KR 
Planck brane TeV brane 



Fig. 7. Warped extra dimension. Even tliougli the standard model particles are shown to be on the 
TeV brane, they may propagate in the bulk depending on the models. 



the Fourier-mode expansion for the mode ipn ix)e ^^v/^ gives 

(*7^5^ - «^7') 4'n{x) = 0. (3.9) 

After a chiral rotation ipn V'ne*^'^ ^^"^n, the second term turns into the usual 
mass term without 7^. The problem is that here are always two eigenvalues 7^ = 
±1 and you find both left- and right-handed fermions with the same quantum 
numbers. Namely, you get Dirac fermion, not Weyl fermion. Then you don't get 
the standard model that distinguishes left from right. In terms of spectrum, what 
is on the left of Fig.[8]is the spectrum because the Fourier modes n and —n give 
the degenerate mass n/R each of them with its own Dirac fermion. 

JL 

4/R n=±4 n=4 

3/R n=±3 n=3 

2/R n=±2 n=2 

1/R n=±1 n=1 

" " ^ ^ , r, n=0 — „ n=0 
L R L R L R 

Fig. 8. The spectrum of fermions in the 5D bulk. After orbifold identification in Fig.|9] the spectrum 
is halved and one can obtain chiral fermions in 4D. 

The trick to get chiral fermions is to use an orbifold Fig.|9] Out of a circle 
(y £ [— Ti'-R, Tri?]), you identify points y and — y to get a half-circle 3^/7^2- There 
are two special points, y = and ttR, that are identified only with themselves 
called "fixed points." In addition, we take the boundary condition that ^{y) — 
— 7^V(~2/)- For n 7^ 0, we use cosny/R for 7^ — —1 and s'mny/R for 7^ = 
+1, without the degeneracy between n and —n. For n = 0, only 7^ = — 1 
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survives with the wave function V'o(?/) = 1- This way, we keep only a half of the 
states as shown in Fig. [8] and we can get chiral fermions. As a consequence, we 
find the system to have a Z2 symmetry under y ttR — y, under which modes 
with even n are even and odd n odd. This Z2 symmetry is called KK parity and 
the lightest KK state (LKP) becomes stable. At the tree-level, all first Kaluza- 
Klein states are degenerate toi = 1/rE Radiative corrections split their masses, 
and typically the first Kaluza-Klein excitation of the U{1)y boson is the LKP 
1461. Because the mass splittings are from the loop diagrams, they are small. 
Similarly to supersymmetry, there is a large number of new particles beyond the 
standard model, namely Kaluza-Klein excitations. Its collider phenomenology 
very much resembles that of supersymmetry and it is not trivial to tell them apart 
at the LHC (dubbed "bosonic supersymmetry" |47|). 




Fig. 9. The orbifold /Z,2. Points connected by tlie solid lines are identified. 



4. Evidence for Dark Matter 



Now we turn our attention to the problem of non-baryonic dark matter in the 
universe. Even though this is a sudden change in the topic, you will see soon that 
it is connected to the discussions we had on physics beyond the standard model. 
We first review basics of observational evidence for non-baryonic dark matter, 
and then discuss how some of the interesting candidates are excluded. It leads 
to a paradigm that dark matter consists of unknown kind of elementary particles. 
By a simple dimensional analysis, we find that a weakly coupled particle at the 
TeV-scale naturally gives the correct abundance in the current universe. We will 
take a look at a simple example quite explicitly so that you can get a good feel 
on how it works. Then I will discuss more attractive dark matter candidates that 
arise from various models of physics beyond the standard model I discussed in 
the previous section. 



'^Here I've ignored possible complications due to brane operators and electroweak symmetry 
breaking. 
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The argument for the existence of "dark matter," namely mass density that 
is not luminous and cannot be seen in telescopes, is actually very old. Zwicky 
back in 1933 already reported the "missing mass" in Coma cluster of galaxies. 
By studying the motion of galaxies in the cluster and using the virial theorem 
(assuming of course that the galactic motion is virialized) he determined the mass 
distribution in the cluster and reported that a substantial fraction of mass is not 
seen. Since then, the case for dark matter has gotten stronger and stronger and 
most of us regard its existence established by now. I refer to a nice review for 
details |48| written back in 1997, but I add some important updates since the 
review. 

Arguably the most important one is the determination of cosmological param- 
eters by the power spectrum of CMB anisotropy. In the fit to the power-law flat 
ACDM model gives flnh^ = 0.127t^;J^?^ and rish^ = 0.0223ti];^!^Sj^ ll49l . 
The point here is that these two numbers are different. Naively subtracting the 
baryon component, and adding the errors by quadrature, I find {flM — ^B)h^ = 
O.lOS^Q g^g 7^ at a very high precision. This data alone says most of the matter 
component in the universe is not atoms, something else. 

Another important way to determine the baryon density of the universe is 
based on Big-Bang Nucleosynthesis (BBN). The baryon density is consistent 
with what is obtained from the CMB power spectrum, fis/i^ = 0.02161q;qo2i 
from five best measurements of deuterium abundance [SO] using hydrogen gas 
at high redshift (and hence believed to be primordial) back-lit by quasars. This 
agrees very well with the CMB result, even though they refer to very different 
epochs: T ~ 1 MeV for BBN while T 0.1 eV for CMB. This agreement gives 
us confidence that we know Vl b very well. 

A novel technique to determine Q,m uses large-scale structure, namely the 
power spectrum in galaxy-galaxy correlation function. As a result of the acoustic 
oscillation in the baryon-photon fluid, the power spectrum also shows the "baryon 
oscillation" which was discovered only the last year |5r|. Without relying on the 
CMB, they could determine VlMh"^ ~ 0.130 ± 0.010. Again this is consistent 
with the CMB data, confirming the need for non-baryonic dark matter 

I'd like to also mention a classic strong evidence for dark matter in galaxies. 
It comes from the study of rotation curves in spiral galaxies. The stars and gas 
rotate around the center of the galaxy. For example, our solar system rotates in 
our Milky Way galaxy at the speed of about 220 km/sec. By using Kepler's law, 
the total mass M(r) within the radius r and the rotation speed at this radius v(r) 
are related by 



v{rf = Gn 



M{r) 



(4.1) 



r 



Once the galaxy runs out of stars beyond a certain r, the rotation speed is hence 
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expected to decrease as v{r) cx r^^/^. This expectation is not supported by 
observation. 

You can study spiral galaxies which happen to be "edge-on." At the outskirts 
of a galaxy, where you don't find any stars, there is cold neutral hydrogen gas. 
It turns out you can measure the rotation speed of this cold gas. A hydrogen 
atom has hyperfine splitting due to the coupling of electron and proton spins, 
which corresponds to the famous A = 21 cm line emission. Even though the 
gas is cold, it is embedded in the thermal bath of cosmic microwave background 
whose temperature 2.7 K is hot compared to the hyperfine excitation hc/kX = 
0.069 kEI Therefore the hydrogen gas is populated in both hyperfine states and 
spontaneously emits photons of wavelength 21 cm by the Ml transition. This can 
be detected by radio telescopes. Because you are looking at the galaxy edge-on, 
the rotation is either away or towards us, causing Doppler shifts in the 21 cm line. 
By measuring the amount of Doppler shifts, you can determine the rotation speed. 
Surprisingly, it was found that the rotation speed stays constant well beyond the 
region where stars cease to exist. 
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Fig. 10. Rotation curve of a spiral galaxy 1521 . 

I mentioned this classic evidence because it really shows galaxies are filled 
with dark matter. This is an important point as we look for signals of dark matter 
in our own galaxy. It is not easy to determine how much dark matter there is, 
however, because eventually the hydrogen gas runs out and we do not know how 
far the flat rotation curve extends. Nonetheless, it shows the galaxy to be made 
up of a nearly spherical "halo" of dark matter in which the disk is embedded. 




'^If we had lived in a universe a hundred times larger, we would have lost this opportunity of 
studying dark matter content of the galaxies ! 
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5. What Dark Matter Is Not 

We don't know what dark matter is, but we have learned quite a bit recently 
what it is not. I have already discussed that it is not ordinary atoms (baryons). 
I mention a few others of the excluded possibilities. 

5.1. MACHOs 

The first candidate for dark matter that comes to mind is some kind of astronom- 
ical objects, namely stars or planets, which are is too dark to be seen. People 
talked about "Jupiters," "brown dwarfs," etc. In some sense, that would be the 
most conservative hypothesisl^ Because dark matter is not made of ordinary 
atoms, such astronomical objects cannot be ordinary stars either. But one can 
still contemplate the possibility that it is some kind of exotic objects, such as 
black holes. Generically, one refers to MACHOs which stand for MAssive Com- 
pact Halo Objects. 

Black holes may be formed by some violent epochs in Big Bang (primordial 
black holes or PBHs) 1531 (see also 1541). If the entire horizon collapses into a 
black hole, which is the biggest mass one can imagine consistent with causality, 
for example in the course of a strongly first order phase transition, the black hole 
mass would be 



Therefore, there is no causal mechanism to produce PBHs much larger than 
IQ^Mq assuming that universe has been a normal radiation dominated universe 
for T < 3 MeV to be compatible with Big-Bang Nucleosynthesis. Curiously, one 
finds A/pBH ~ Mq if it formed at the QCD phase ti-ansition T w 100 MeV f55l. 
On the other hand, PBHs cannot be too small because otherwise they emit Hawk- 
ing radiation of temperature T = (87rGArAfpBH)^ that would be visible. The 
limit from diffuse gamma ray background implies A/pbh ^ 1O~^^M0. 

How do we look for such invisible objects? Interestingly, it is not impossible 
using the gravitational microlensing effects ll56ll . The idea is simple. You keep 
monitoring millions of stars in nearby satellite galaxies such as Large Magel- 
lanic Cloud (LMC). Meanwhile MACHOs are zooming around in the halo of our 
galaxy at u « 220 km/s. By pure chance, one of them may pass very close along 
the line of sight towards one of the stars you are monitoring. Then the gravity 
would focus light around the MACHO, effectively making the MACHO a lens. 
You typically don't have a resolution to observe distortion of the image or multi- 
ple images, but the focusing of light makes the star appear temporarily brighter. 

"Somehow I can't call primordial black holes a "conservative" candidate without chuckling. 




(5.1) 
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This is called "microlensing." By looking for such microlensing events, you can 
infer the amount of MACHOs in our galactic halo. 

I've shown calculations on the deflection angle by the gravitational lensing 
and the amplification in the brightness in the appendix. (Just for fun, I've also 
added some discussions on the strong lensing effects.) The bottom line is that 
you may expect the microlensing event at the rate of 



rate « 5 x 10 



1/2 



year V^macho/ 
towards the LMC, with the duration of 



(5.2) 



duration ; 



6 X lO'^sec ( MM^) 
\ Me J 



1/2 



25kpc 



(5.3) 



where c?i (1^2) is the distance between the MACHO and us (the lensed star). 

Two collaborations, the MACHO collaboration and the EROS collaboration, 
have looked for microlensing events. The basic conclusion is that MACHOs of 
mass 10^^-30 Mq cannot make up 100% of our galactic halo (Fig.fTTTi. See also 
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Fig. 11. Limit on the halo fraction / of MACHOs from the EROS collaboration 1571 . The spherical 
isothermal model of halo predicts the optical depth towards the LMC of t = 4.7 X 10"'^. For more 
details, see the paper 



Even though the possibility of MACHO dark matter may not be completely 
closed, it now appears quite unlikely. The main paradigm for the dark matter of 
the universe has shifted from MACHOs to WIMPs. 
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5.2. Neutrinos 

Having discovered neutrinos have finite mass, it is also natural to consider neu- 
trinos to be dark matter candidate. As a matter of fact, neutrinos are a component 
of dark matter, contributing 



~ 94eV ■ ^ ' 



It is an attractive possibility if the particles which we already know to exist could 
serve as the required non-baryonic dark matter. 

However, as Sergio Pastor discussed in his lectures, neutrinos are not good 
candidates for the bulk of dark matter for several reasons. First, there is an upper 
limit on neutrino mass from laboratory experiments (tritium beta decay) m < 
2 eV ||59ll . Combined with the smallness of mass-squared differences Attiq = 
8 X 10~^ eV^ and Am^ = 2.5 x 10"'^ eV^, electron-volt scale neutrinos should 
be nearly degenerate. Then the maximum contribution to the matter density is 
n^h'^ < (3 X 2/94) < 0.064. This is not enough. 

Second, even if the laboratory upper limit on the neutrino mass turned out to be 
not correct, there is a famous Tremaine-Gunn argument (60\. For the neutrinos 
to dominate the halo of dwarf galaxies, you need to pack them so much that 
you would violate Pauli exclusion principle. To avoid this, you need to make 
neutrinos quite massive > 500 eV so that you don't need so many of them ||6T1 . 
This obviously contradicts the requirement that il^ < 1. 

Third, neutrinos are so light that they are still moving at speed of light (Hot 
Dark Matter) at the time when the structure started to form, and erase structure at 
small scales. Detailed study of large scale structure shows such a hot component 
of dark matter must be quite limited. The precise limit depends on the exact 
method of analyses. A relatively conservative limit says J^i '^I'i < 0-62 eV 1*621 
while a more aggressive limit goes down to 0.17 eV |63|. Either way, neutrinos 
cannot saturate what is needed for non-baryonic dark matter 

In fact, what we want is Cold Dark Matter, which is already non-relativistic 
and slowly moving at the time of matter-radiation equality T ^ 1 eV. Naively a 
light (sub-electronvolt) particle would not fit the bill. 

A less conservative hypothesis may be to postulate that there is a new heavy 
neutrino (4th generation). This is a prototype for WIMPs that will be discussed 
later. It turns out, however, that the direct detection experiments and the abun- 
dance do not have a compatible mass range. Namely the neutrinos are too strongly 
coupled to be the dark matter! 
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5.3. CHAMP s and SIMP s 

Even though people do not talk about it any more, it is worth recalling that dark 
matter is unlikely be charged (CHAMP) |64| or sti-ongly interacting (SIMP) l!65l . 
I simply refer to papers that limit such possibilities, from a multitude of search 
methods that include search for anomalously heavy "water" molecule in the sea 
water, high-energy neutrinos from the center of the Earth from annihilated SIMPs 
accumulated there, collapsing neutron stars that accumulate CHAMPs, etc0 

6. WIMP Dark Matter 

WIMP, or Weakly Interactive Massive Particle, is the main current paradigm for 
explaining dark matter of the universe. With MACHOs pretty much gone, it is 
indeed attractive to make a complete shift from astronomical objects as heavy as 
Mq « 10^^ GeV to "heavy" elementary particles of mass ^ 10^ GeV. I will 
discuss why this mass scale is particularly interesting. 

6.1. WIMP 

The idea of WIMP is very simple. It is a relatively heavy elementary particle x 
so that accelerator experiments so far did not have enough energy to create them, 
namely ^10^ GeV. On the other hand, the Big Bang did once have enough 
energy to make them. 

Let us follow the history from when T > m^. WIMPs were created as much 
as any other particles. Once the temperature dropped below m^, even the uni- 
verse stopped creating them. If they are stable, whatever amount that was pro- 
duced was there, and the only way to get rid of them was to get them annihilating 
each other into more mundane particles {e.g., quarks, leptons, gauge bosons). 
However, the universe expanded and there were fewer and fewer WIMPs in a 
given volume, and at some point WIMPs stopped finding each other. Then they 
could not annihilate any more and hence their numbers become fixed ("freeze- 
out"). This way, the universe could still be left with a certain abundance of 
WIMPs. This mechanism of getting dark matter is called "thermal relics." 

Let us make a simple estimate of the WIMP abundance. In radiation domi- 
nated universe, the expansion rate is given by 



I once got interested in tlie possibility that Jupiter is radiating lieat more tlian it receives from 
the Sun because SIMPs are annihilating at its core |66|. It does not seem to explain heat from other 
Jovian planets, however, once empirical limits on SIMPs are taken into account. 




(6.1) 
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where Mpi = l/^/EnU^ = 2.4 x 10^^ GeV is the reduced Planck scale. For 
simple estimates, we regard (7r^/90)^/^ = 0.33 ~ 1 and ignore many other fac- 
tors of 0(1). Hence, H ~ gV^T"^ /Mpi. The entropy density is correspondingly 



27r2 



2\ 1/2 



s = g,T^ y—j ~ g,T\ (6.2) 

Given the thermally averaged annihilation cross section {aannv), and the number 
density of WIMPs ri;^, the annihilation rate of a WIMP is 

r = {(TannV)n^. (6.3) 

The annihilation stops at the "freeze-out temperature" Tf when F ~ iJ, and 
hence 

n^{Tf)^g'J' ^ . (6.4) 

The yield of WIMPs is defined by = n^/s. This is a convenient quantity 
because it is conserved by the expansion of the universe as long as the expansion 
is adiabatic, i.e., no new source of heat. This is due to the conservation of both 
the total entropy and total number of particles and their densities both scale as 
The estimate of the yield is 



V -1/2 1 -1/2 

Yy g^ = g'* . (6.5) 

{crannV)TfMpi {cTannV)m^Mpi 



Here, we defined Tj = m^/xf. We will see later from more detailed calculations 
that X/ ~ 20. The abundance in the current universe is calculated using the yield 
and the current entropy density, divided by the current critical density. 

So -1/2 Xf So 

S Pc {(TannV)Mpi 

We use So = 2890 cm^'^ and pc — 1.05 x 10~^/i^ GeV cm~^, where the current 
Hubble constant is Hq = lOO/i km/sec/Mpc with h « 0.65. In order of obtain 
f2;^/i^ ~ 0.12, we need 

. ,:"^,,hll21^p^ . 10- GeV-. ,6.7, 

Recall a typical annihilation cross section of a particle of mass by a relatively 
weak interaction of electromagnetic strength (e.g., e+e~ — > 77) is 

ml 
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To obtain the correct abundance, what we need is 

~ 300 GeV. (6.9) 

This is a very interesting result. Namely, the correct abundance of thermal relics 
is obtained for a particle mass just beyond the past accelerator Umits and where 
we expect new particles to exist because of the considerations of electroweak 
symmetry breaking and the hierarchy problem. In other words, it is exactly the 
right mass scale for a new particle! 

In the next few sections, we will firm up this naive estimate by solving the 
Boltzmann equations numerically. We will also study a concrete model of a new 
particle for dark matter candidate and work out its annihilation cross section. In 
addition, we will see if we have a chance of "seeing" the dark matter particle 
in our galactic halo, or making it in future accelerator experiments. Then we 
will generalize the discussions to more theoretically attractive models of physics 
beyond the standard model. 



6.2. Boltzmann Equation 

You have already seen Boltzmann equation in lectures by Sabino Matarrese and 
I don't repeat its derivations. We assume kinetic equilibrium, namely that each 
particle species has the Boltzmann distribution in the momentum space except 
for the overall normaUzation that is given by its number density. Considering the 
process of X1X2 ^ XsXi^ where Xi refers to a certain elementary particle, the 
Boltzmann equation for the number density ni for the particle xi is 

Here, av is the cross section common for the process X1X2 X3X4 and its in- 
verse process XsXi ~^ X1X2 assuming the time reversal invariance. The number 
densities with superscript ° refer to those in the thermal equiUbrium. 

In the case of our interest, xsA are "mundane" Ught (relativistic) particles in 
the thermal bath, and hence 713,4 = 713 4. In addition, we consider the annihilation 
XX (mundane)^, and hence ni = 712. The Boltzmann equation simplifies 
drastically to 

a-'^ = Knnv)[{n'^^r - {n^n (6.11) 



This time we pay careful attention to all numerical factors. We use 
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9*T' r-^] , (6.13) 



45 



87r„ 7r2„, ^2 y4 

72 

X ^ (6.15) 

Even though we start out at temperatures T > when x relativistic, 
eventually the temperature drops below and we can use non-relativistic ap- 
proximations. Then the equilibrium number density can be worked out easily 
as 

,0 ^ /^e-^/T (E = rn ' ^ 



27^ J ^ " {2Trxf/ 



= e-^^-{!I^X'' = e-^,^. (6.16) 



Therefore 



Changing the variables from to Y and t to the Boltzmann equation 
becomes 

Here, we used s(r) = s{m^)/x^ and 

1 m2 1 
dt = —F7rf\^ = - rrf \rrz '^T=— -xdx. (6.19) 

It is useful to work out 

H{m^) 45 \tt 



^ gj m-^Mpi = 1.32g:'m^Mpi. (6.20) 



Note that the annihilation cross section (cTannw) is insensitive to the temperature 
once the particle is non-relativistic T <C m^l^ Therefore the whole combina- 
tion (o'annv) is just a dimensionless number. The only complication is 



''This statement assumes that the annihilation is in the S-wave. If it is in the /-wave, {cFannv) on 
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that Yq has a strong dependence on x. We can further simphfy the equation by 
introducing the quantity 



We obtain 

i^ = --,{v'-yl). (6.22) 
ax x'^ 

with 

yo = 0.192gr'/'Mpim;,(a,„„?;)a:3/2e-^ (6.23) 



6.5. Analytic Approximation 

Here is a simple analytic approximation to solve Eq. ( 16.221 ). We assume Y tracks 
Yq for X < Xf. On the other hand, we assume y ^ yo for x > xj because yo 
drops exponentially as e^^. Of course this approximation has a discontinuity at 
X = Xf, but the transition between these two extreme assumptions is so quick 
that it turns out to be a reasonable approximation. Then we can analytically solve 
the equation for x > xj and we find 

' ' ' (6^24, 



y{oo) y{xf) Xf ' 

Since y(oo) ^ y{xf), we obtain the simple estimate 

y{(X))^Xf. (6.25) 

Given this result, we can estimate Xf as the point where yo {x) drops down ap- 
proximately to Xf, 

0.192g;^^h4pim^{aannv)xf\-''f « Xf, (6.26) 



and hence 

Xf « 



In 



/ 1 /2 

0.192m^Mpi {aannV)Xj:' 



1/2 

9* 



24 + In + In ^"""""^ , - 1 In ^. (6.27) 

100 GeV 10-9GeV-2 2 100 
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6.4. Numerical Integration 

I've gone through numerical integration of the Boltzmann equation Eq. (16.221 1. 
Fig. [12] shows the x-evolution of y. You can see that it traces the equilibrium 
value very well early on, but after x of about 20, it starts to deviate significantly 
and eventually asymptotes to a constant. This is exactly the behavior we expected 
in the analytic approximation studied in the previous section. 

10^0 

y 10^ 

10'' 
100 



10 100 1000 10000 
X 

Fig. 12. Numerical solution to the Boltzmann equation Eq. 16. 22) for m = 100 GeV, = 100, 
(f ann'f) = 10~® GeV~^. Superimposed is the equilibrium value yo. 

Fig.[T3]shows the asymptotic values y(oo) which we call Xf. I understand this 
is a confusing notation, but we have to define the "freeze-out" in some way, and 
the analytic estimate in the previous section suggest that the asymptotic value 
y{oo) is nothing but the freeze-out value Xf. This is the result that enters the 
final estimate of the abundance and is hence the only number we need in the 
end anyway. It does not exactly agree with the estimate in the previous section, 
but does very well once I changed the offset in Eq. (I6.27l i from 24 to 20.43. 
Logarithmic dependence on is verified beautifully. 

24 
23 
22 
21 
20 

100 200 500 1000 2000 

(GeV) 

Fig. 13. Xf values as afunction of m;,;., for 3, = 100, (uannv) = 10~^ GeV~^. The dots are the 
results of numerical integrations, while the solid line is just In m-^ with an offset so that Xf = 20.43 
for my = 100 GeV. 
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Flitting everything back together, 

Him^) Xf 

Py = m^ny = m^Ys = — — -s (6.28) 

'^^ ^ ^ ^ X (cTannV) 

As before, we use sq — 2890 cm^^ and pc — 1.05 x 10~^/i^ GeV cm^"^, where 
the current Hubble constant is Hq = IQQh km/sec/Mpc with h « 0.65. To obtain 
^Mh^ — 0.12, we find (cra„„f) = 1.6 x 10^^ GeV^^, confirming the simple 
estimate in Section|6Tl 



6.5. The New Minimal Standard Model 

Now we would like to apply our calculations to a specific model, called the New 
Minimal Standard Model |2 |. This is the model that can account for the empirical 
facts listed in Section 12.11 with the minimal particle content if you do not pay 
any attention to the theoretical issues mentioned in Section |272l It accomplishes 
this by adding only four new particles to the standard modelllj very minimal 
indeed! The dark matter in this model is a real scalar field S with an odd Z2 
parity S ^ —S, and its most general renormalizable Lagrangian that should be 
added to the Standard Model Eq. ( 12. Il l is 

Cs = \d,Sd^S - \mlS' - \\H\'S^ - (6.29) 

The scalar field S is the only field odd under Z2, and hence the S boson is stable. 
Because of the analysis in the previous sections, we know that if ms is at the 
electroweak scale, it may be a viable dark matter candidate as a thermal relic. 
This is a model with only three parameters, ms, k, and h, and actually the last 
one is not relevant to the study of dark matter phenomenology. Therefore this is 
a very predictive model where one can work it out very explicitly and easily. 

To calculate the dark matter abundance, what we need to know is the annihila- 
tion cross section of the scalar boson 5*. This was studied first in |68 1 and later in 
||69l , but the third diagram was missing. In addition, there is a theoretical bounds 
on the size of couplings k and h so that they would stay perturbative up to high 
scales {e.g., Planck scale). The cosmic abundance is determined by ms and k in 
addition to m^- Therefore on the {k^mh) plane, the correct cosmic abundance 
determines what ms should be. This is shown in Fig. [15] You can see that for a 
very wide range 7715 ~ 5.5 GeV-1.8 TeV, the correct cosmic abundance can be 
obtained within the theoretically allowed parameter space. For heavy ^ m^, 
the cross section goes like k"^ /m^g and is independent of mh. This is why the ms 
contours are approximately straight vertically. For light ms ^ mu, the cross 



'^The other three are the inflaton and two right-handed neutrinos. 
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S f S V = W,Z 




S h S h s 



Fig. 14. Feynman diagrams for the anniliilation of S scalars. Tlie final states in the first diagram can 
be any of the quark or lepton pairs //. 




) h . i . i:. I.-. . , , I , , . / I . , , . I , . , . li 
0.1 0.2 0.3 0.4 0.5 



Fig. 15. The region of the NMSM parameter space (A;(mz),r?i(j) that satisfies the stability and 
triviality bounds, for h(mz) = 0, 1.0, and 1.2. Also the preferred values from the cosmic abundance 
Qsh^ = 0.11 are shown for various mg. Taken from (2). 



section goes like k'^rn^/m'^. This is why the ms contours approximately have 
a fixed k/m\ ratio. Note that when ms — m/i/2, the first two diagrams can hit 
the Higgs pole and the cross section can be very big even for small k. This res- 
onance effect is seen in Fig. [T5] where ms = 75 GeV line reaches almost A: ~ 
form,, = 150 GeV. 

You may wonder why I am talking about S as light as 5.5 GeV. Shouldn't 
we have seen it already in accelerator experiments? Actually, no. The only 
interaction the S boson has is with the Higgs boson which we are yet to see. 
Therefore, we could not have produced the S boson unless we had produced the 
Higgs boson. That is why even such a light 5* boson does not contradict data. In 
other words, it wouldn't be easy to find this particle in accelerator experiments. 
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6.6. Direct Detection Experiments 

How do we know if dark matter is indeed in the form of WIMP candidate you 
like? One thing we'd love to see is the direct detection of WIMPs. The idea is 
very simple. You place a very sensitive device in a quite location. WIMPs are 
supposed to be flying around in the halo of our galaxy with the typical speed of 
~ 220 km/s ^ Because they are only very weakly interacting, they can go 

through walls, rocks, even the entire Earth with little trouble, just like neutrinos. 
For a mass of ^ 100 GeV, its typical kinetic energy is i?kin = ^"T'x^^ ~ 
50 keV. If the WIMP (ever) scatters off an atomic nucleus, the energy deposit is 
only (at most) of this order of magnitude. It is a tiny energy deposit that is very 
difficult to pick out against background from natural radioactivity (typically Me V 
energies). Therefore you have to make the device very clean, and also place it 
deep underground to be shielded from the cosmic -ray induced backgrounds, most 
importantly neutrons ejected from the rocks by cosmic -ray muons. One you've 
done all this, what you do is to wait to see this little "kick" in your detector 

Let us do an order of magnitude estimate. The local halo density is esti- 
mated to be about w 0.3 GeV/cm^. The number density of WIMPs is 
Uy, — p^^^°/m^. The flux of WIMPs is roughly vuy, « 10~^cn^. The elastic 
cross section of WIMP on neutron or proton may be spin-independent or spin- 
dependent. In the spin-independent case, the amplitude of the WIMP-nucleus 
cross section goes as A (mass number) and hence the cross section on the nucleus 
a A goes as A^. Of course the detailed scaling is model-dependent, but in most 
phenomenological analyses (and also analyses of data) we assume aA = A^ap. 
Let us also assume ^^Ge as the detector material so that A = 56. Then the 
expected event rate is 

^ ^halo^. ^'ttarget ^^ _ 10 100 GcV mtarget ^ CTp 

^ ruA year 100 kg 56 lO^'*^ cm^ 

To prepare a very sensitive device as big as 100 kg and make it very clean is a 
big job. You can see that your wait may be long. 

Now back to the New Minimal Standard Model. The scattering of the 5* boson 
off a proton comes from the t-channel Higgs boson exchange. The coupling of 
the Higgs boson to the nucleon is estimated by the famous argument |70| using 
the conformal anomaly. The mass of the proton is proportional to the QCD scale 
rup cx Aqcd (m-u, rrid, rris are ignored and hence this is the three-flavor scale). 
It is related to the Higgs expectation value through the one-loop renormalization 
group equation as (we do not consider higher loop effects here) 

Aqcd = m^/^m^/^^m^/^M^e-s^'/^'^^^) (6.31) 
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where M is some high scale and each quark mass is proportional to v. The 
coupling of the Higgs to the proton is given by expanding the vacuum expectation 
value as V V + h, and hence 

dmp 2 TOp 

^-'^ = ^ = 9 V ^'-^'^ 

This allows us to compute the scattering cross section of the S boson and the 
nucleon. 




Fig. 16. Feynman diagrams for the scattering of the S-boson off a proton or neutron. 



The result is shown in Fig.[T6]as the red solid line. The point is that the elastic 
scattering cross section tends to be very small. Note that a hypothetical neutrino 
of the similar mass would have a cross section of dp ~ Gj^ml/n ~ lO"^** cm^ 
which is much much bigger than this. This is the typical WIMP scattering cross 
section. Superimposed is the limit from the CDMS-II experiment ITTTI and hence 
the direct detection experiments are just about to reach the required sensitivity. In 
other words, this simple model is completely viable, and may be tested by future 
experiments. For the resonance region nis — rah/2, the coupling k is very small 
to keep enough abundance and hence the direct detection is very difficult. 

The future of this field is not only to detect WIMPs but also understand its 
identity. For this purpose, you want to combine the accelerator data and the 
direct detection experiments. The direct detection experiments can measure the 
energy deposit and hence the mass of the WIMP. It also measures the scattering 
cross section, even though it suffers from the astrophysical uncertainty in the 
estimate of the local halo density. On the other hand, assuming ms < mh/2, 
the Higgs boson decays invisibly h — > SS. Such an invisible decay of the light 
Higgs boson can be looked for at the LHC using the Vl^-fusion process. Quarks 
from both sides radiate an off-shell -boson that "fuse" in the middle to produce 
a Higgs boson. Because of the kick by the off-shell Ty-boson, the quarks acquire 
Pt ~ m-w l'2 in the final state and can be tagged as "forward jets." Even though 
the Higgs boson would not be seen, you may "discover" it by the forward jets 
and missing Et |74|. The ILC can measure the mass of the Higgs precisely 
even when it decays dominantly invisibly (see, e.g., 1751 ) and possibly its width. 
Combining it with the mass from the direct detection experiments, you can infer 
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Fig. 17. The elastic scattering cross section of Dark Matter from nucleons in NMSM, as a function 
of the Dark Matter particle mass mg for rrtj, = 150 GeV. Note that the region mg > 1.8 TeV is 
disallowed by the triviality bound on k. Also shown are the experimental bounds from CDMS-II 1711 
and DAMA [72], as well as improved sensitivities expected in the future |73 1. Taken from (2). 



the coupling k and calculate its cosmic abundance. It would be a very interesting 
test if it agrees with the cosmological data Q,m ~ 0.23. If it does, we can claim 
a victory; we finally understand what dark matter is\ 

6.7. Popular WIMPs 

Superpartners of the photon and Z, and neutral Higgs bosons (there are two of 
them), mix among each other once SU{2) x U (1) symmetry breaks. Out of four 
such "neutralino" states, the lightest one is often the LSlO and is the most pop- 
ular candidate for dark matter in the literature (see, e.g., IITTl l78l |79 1 for some of 
the recent papers). One serious problem with the supersymmetric dark matter is 
that there are many parameters in the model. Even the Minimal Supersymmet- 
ric Standard Model (MSSM) has 105 more parameters than the standard model. 
It is believed that the fundamental theory determines all these parameters or at 
least reduces the number drastically, and one typically ends up with five or so 
parameters in the study. Depending on what parameter set you pick, the phe- 
nomenology may be drastically different. For the popular parameter set called 
CMSSM (Constrained MSSM, also called minimal supergravity or mSUGRA) 
with four parameters and one sign, see a recent study in |76|. I do not go into 



Superpartner of neutrinos is not out of the question |80| . 
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detailed discussions about any of them. I rather mention a few generic points. 

First, we do get viable dark matter candidates from sub-TeV supersymmetry 
as desired by the hierarchy problem. This is an important point that shouldn't 
be forgotten. Second, what exactly is the mass and composition of the neu- 
tralino depends on details of the parameter set. The supersymmetric standard 
model may not be minimal either; an extension with additional singlet called 
Next-to-Minimal Supersymmetric Standard Model (NMSSM) is also quite pop- 
ular Third, sub-TeV supersymmetry can be studied in great detail at LHC and 
(hopefully) ILC, so that we can measure their parameters very precisely (see [811 
for the early work). One can hope to correlate the accelerator and underground 
data to fully test the nature of the dark matter |82|. Fourth, a large number of 
particles present in the model may lead to interesting effects we did not consider 
in the discussions above. I'll discuss one of such effects briefly below. 

Universal Extra Dimension (UED) is also a popular model. I'm sure Geraldine 
Servant will discuss dark matter in this model in her lectures because she is one 
of the pioneers in this area ll86l . Because the LKP is stable, it is a dark matter 
candidate. Typically the first KK excitation of the U{1)y gauge boson is the 
LKP, and its abundance can be reasonable. Its prospect for direct and indirect 
detection experiments is also interesting. See a review article that came out after 
Les Houches fWi\. 

The most striking effect of having many particle species is the coannihilation 
ll83l . An important example in the case of supersymmetry is bino and stau. Brno 
is the superpartner of the U{\)y gaugino (mixture of photino and zino), and its 
annihilation cross section tends to be rather small partly because it dominantly 
goes through the P-wave annihilation. If, however, the mass of the stau is not too 
far above bino, stau is present with the abundance suppressed only by e~^™/^ 
(Am — iTif — nig) assuming they are in chemical equilibrium Bt ^f. There 
are models that suggest the mass splitting is indeed small. The cross sections 
Bt — > 7T and ff* 77, // etc tend to be much larger, the former going 
through the S'-wave and the latter with many final states. Therefore, despite the 
Boltzmann suppression, these additional contributions may win over (a-ggv). 
There are other cases where the mass splitting is expected to be small, such as 
higgsino-like neutralinos [84]. In the UED, the LKP is quite close in its mass to 
the low-lying KK states and again coannihilation is important. 

6.8. Indirect Detection Experiments 

On the experimental side, there are other possible ways of detecting signals of 
dark matter beyond the underground direct detection experiments and collider 
searches. They are indirect detection experiments, namely that they try to detect 
annihilation products of dark matter, not the dark matter itself. For annihilation to 
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occur, you need some level of accumulation of dark matter The possible sites are: 
galactic center, galactic halo, and center of the Sun. The annihilation products 
that can be searched for include gamma rays from the galactic center or halo, 
e+ from the galactic halo, radio from the galactic center, anti-protons from the 
galactic halo, and neutrinos from the center of the Sun. Especially the neutrino 
signal complements the direct detection experiments because the sensitivity of 
the direct detection experiments goes down as 1 /m^ because the number density 
goes down, while the sensitivity of the neutrino signal remain more or less flat 
for heavy WIMPs because the neutrino cross section rises as cx m^. You can 
look at a recent review article lf85l on indirect searches. 



7. Dark Horse Candidates 



7.1. Gravitino 



Assuming i?-parity conservation, superparticles decay all the way down to what- 
ever is the lightest with odd i?-parity. We mentioned neutralino above, but an- 
other interesting possibility is that the LSP is the superpartner of the gravitino, 
namely gravitino. Since gravitino couples only gravitationally to other particles, 
its interaction is suppressed by 1/Mpi. It practically removes the hope of direct 
detection. On the other hand, it is a possibility we have to take seriously. This is 
especially so in models with gauge-mediated supersymmetry breaking 1881 . 
The abundance of light gravitinos is given by 

^3/2^^ = — ^ (7.1) 
2keV ^ ' 

if the gravitinos were thermalized. In most models, however, the gravitino is 
heavier and we cannot allow thermal abundance. The peculiar thing about a 
light gravitino is that the longitudinal (helicity ±1/2) components have a much 
stronger interaction TO5c/5y/(TO3/2-/W^Pi) if "^3/2 *C rnsuSY- Therefore the 
production cross section scales as ct ~ ^"susy/ ^^'^ 'he abundance 
scales as Y3/2 - aT^/H{T) - mljjgyT / [mlf^Mpi). Therefore we obtain 
an upper limit on the reheating temperature after the inflation ll89l l90l . If the 
reheating temperature is right at the limit, gravitino may be dark matter. 

There is however another mechanism of gravitino production. The abundance 
of the "LSP"is determined the usual way as a WIMP, while it eventually decays 
into the gravitino. This decay can upset the Big-Bang Nucleosynthesis fSOl, but 
it may actually be helpful for a region of parameter space to ease some tension 
among various light element abundances |91 1. Note that the "LSP" (or more cor- 
rectly NLSP: Next-Lightest Supersymmetric Particle) may be even electrically 
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10"'^ 10"^ 10"* 10'' 10'^ 10"' 1 
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Fig. 18. The upper bound on Tmax as a function of the gravitino mass from the requirement that 
the relic stable gravitinos do not overdose the Universe. Taken from |90|. It assumes h = I and 
^3/2^^ < 1 and hence the actual constraint is nearly an order of magnitude more stringent than this 
plot. 



charged or strongly coupled as it is not dark matter. When superparticles are 
produced at colliders, they decay to the "LSP" inside the detector, which escapes 
and most Ukely decays into the gravitino outside the detector. If the "LSP" is 
charged, it would leave a charged track with anomalously high dE/dx. It is in 
principle possible to collect the NLSP and watch them decay, and one may even 
be able to confirm the spin 3/2 nature of the gravitino ll92l and its gravitational 
coupling to matter 193)1 . 

If the gravitino is heavier than the LSP, its lifetime is calcuated as 

r(G ^ 77) = 3.9 x 10^ sec . (7.2) 

It tends to decay after the BBN and upsets its success. Its production cross section 
scales as cT ~ l/Mp^ and hence 1^3/2 ~ T/Mpi. Depending on its mass and 
decay modes, one can again obtain upper limits on the reheating temperature. 
The case of hadronic decay for 7713/2 ~ 1 TeV is the most limiting case that 
requires Tjih ^ 10^ GeV ||94l , causing trouble to many baryogenesis models. 

7.2. Axion 

One of the puzzles about the standard model I discussed earlier is why 6 < 10^^*^ 
in QCD. A very attractive solution to this problem is to promote 6* to a dynamical 
field, so that when it settles to the minimum of the potential, it automatically 
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makes 9 effectively zero 119511 . The dynamical field is called axion, and couples 
(after integrating out some heavy fields) as 



C 



j e'''''"^Tr(G^^Gpcr). 

J a , 



(7.3) 



Here, fa is called axion decay constant which has a dimension of energy. The 
potential is given approximately as 



V 



2 r2 



a 

fa 



(7.4) 



and indeed the axion field settles to a = ^Ofa and the GG term vanishes. The 
axion mass is therefore 



6 /ieV 



10^2 GeV 

fa 



(7.5) 



Various astrophysical limits basically require fa ^ 10^" GeV and hence the 
axion is a very light boson (see, e.g., [96]). Most of these limits come from the 
fact that the axion can carry away energy from stars and would cool them too 
quickly, such as white dwarfs, red giants, and SN1987A. Models of such high 
fa are called "invisible axion" models because then the axion coupling is very 
weak to other particles, avoids these limits, and hence is very difficult to observe. 
There are two popular versions, KSVZ ||97l|98l and DFSZ ISIIOOl models. 

In the early universe T ^ GeV, the axion potential looks so flat that it cannot 
tell where the minimum is E Therefore we expect it starts out wherever it finds 
itself, mostly likely not at the minimum. The likely initial misplacement is of the 
order of fa- Now we would like to know what happens afterwards. 

Let us consider a scalar field in an expanding universe. Neglecting the spatial 
variation and considering the time dependence alone, the equation of motion is 



4> + 3H4) + V'{(l)) =0. 



(7.6) 



For a quadratic potential (mass term) V{(j)) — ^ 
larly simple and homogeneous. 



^m?4P', the equation is particu- 



(7.7) 



'^'This is especially trae for the axion because its mass oiiginates from the QCD instaton effects 
which are suppressed by powers of the temperature in hot thermal bath |101| . 
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It is useful to solve this equation for a constant (time-independent) H first. For 
Fourier modes (f) ^ e~*'^* (of course we take the real part later on), we need to 
solve 

-uj"^ - UH + TO^ = 0, (7.8) 
and we find 

-3ii? ± V-9i?2 + 4^2 



1 



-MH, {H » m) 

m-i^H, -m-i^H {H <^m) 



(7.9) 



Therefore for H ^ m (early universe), one of the solutions damps quickly (j)^ ^ 
g-3Ht^ while the other is nearly stationary ~ e"^'" /3//)t because 
it is "stuck" by the friction term —3H^. On the other hand the field oscillates 
around the minimum as (/> = 0oe~""*e~3-'^*/2. 

One can improve this analysis for time-dependent H{t) ~ l/2t and a{t) oc 
t^^^ in the radiation-dominated univese, by replacing e^*'^±* by e"'-'^*'^*^* 
(adiabatic approximation). When H ^ m {mt <C 1), we find 



t J \ a / 

(p- = (poe •'*o ' = (^oe ' w 00- (7.10) 

The second one is the solution that is stuck by the friction. On the other hand 
when H <^m, 

The field damps as a^'^/^, and its energy as V = 171^(1^/2 cx a""^, just like 
non-relativistic matter. In fact, a coherently oscillating homogeneous field can be 
regarded as a Bose-Einstein condensate of the boson at zero momentum state. 

Therefore, the axion field can sit on the potential and does not roll down be- 
cause of the large friction term —3H(j) when H ^ m. On the other hand for later 
universe H ^ m,\X oscillates as a usual harmonic oscillator e^""* and dilutes 
as non-relativistic matter This is why a very light scalar field can be a candi- 
date for cold dark matter. Counterintuitive, but true. This way of producing cold 
dark matter is called "misaUgnment production" because it is due to the initial 
misalignment of the axion field relative to the potential minimum. Because the 
amount of misalignment is not known, we cannot predict the abundance of axion 
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precisely. Assuming the misalignment of 0{fa), fa — lO^'^ GeV is the preferred 
range for axion dark matter. 

There is a serious search going on for axion dark matter in the halo of our 
galaxy. In addition to the required coupling of the axion to gluons, most models 
predict its coupling to photons a{E ■ B) of the similar order of magnitude. The 
ADMX experiment places a high-Q cavity in a magnetic field. When an axion 
enters the cavity, this coupling would allow the axion to convert to a photon, 
which is captured resonantly by the cavity with a high sensitivity. By changing 
the resonant frequencies in steps, one can "scan" a range of axion mass. Their 
limit has just reached the KSVZ axion model QoSlIIMini and an upgrade to 
reach the DFSZ axion model using SQUID is in the works. See Ml 0511 for more 
on axion microwave cavity searches. 

7.3. Other Candidates 

I focused on the thermal reclic of WIMPs primarily because there is an attractive 
coincidence between the size of annihilation cross section we need for the correct 
abundance and the energy scale where we expect to see new partices from the 
points of view of electroweak symmetry breaking and hierarchy problem. Axion 
is not connected to any other known energy scale, yet it is well motivated from 
the strong CP problem. On the other hand, nature may not necessarily tell us 
a "motivation" for a particle she uses. Indeed, people have talked about many 
other possible candidates for dark matter. You may want to look up a couple of 
keywords: sterile neutrinos, axinos, warm dark matter, mixed dark matter, cold 
and fuzzy dark matter, Q-balls, WIMPZILLAs, etc. Overall, the candidates in 
this list range in their masses from 10^^^ eV to 10^^ eV, not to mention still 
possible MACHOs < lO^^Af© = 10^^ eV. Clearly, we are making progress. 



8. Cosmic Coincidence 

Whenever I think about what the univese is made of, including baryons, photons, 
neutrinos, dark matter, and dark energy, what bothers me (and many other people) 
is this question: why do they have energy densities within only a few orders of 
magnitude? They could have been many orders of magnitude different, but they 
aren't. This question is related to the famous "Why now?" problem. The problem 
is clear in Fig. [191 As we think about evolution of various energy densities over 
many decades of temperatures, why do we live at this special moment when the 
dark matter and dark energy components become almost exactly the same? I feel 



''They ignored theoretical uncertainty in the prediction of the axion-to-photon coupling, and the 
KSVZ model i.s not quite excluded yet [102]. 
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like I'm back to Ptolemy from Copernicus. We are special, not in space any 
more, but in time. Is that really so? 



now 




lO'' 10* 10^ 10° 10^ 10^ 10^ 10^10"™10"^^10""10"'''10""' 



T [GeV] 

Fig. 19. The evolution of radiation, matter, and cosmological constant (A) components of the univese 
as the temperature drops over many orders of magnitude, "Now" is a very special moment when 
matter and A are almost exactly the same, and the radiation is not that different either. 



In Fig. [19] the radiation component goes down as T~'^, while the matter T~^. 
The cosmological constant is by definition constant T". Matter and A meet now. 
When thinking about this problem, it is always tempting but dangerous to bring 
"us" into the discussion. Then we will be forced to talk about conditions for 
emergence of intelligent lifeforms, which we don't know very well about. In- 
stead, it may be better to focus on physical quantities; namely the triple coinci- 
dence problem that three lines with different slopes seem to more or less meet at 
a point. In fact, dimensional analysis based on TeV-scale WIMP suggests 

(Te\^\^ 3 

Pinattcr ^ -n— (8.1) 

which agrees with Pradiation ~ T'^ at the temperature T ~ TeV^/Mp; « meV = 
10 K; this is about now! In order for the cosmological constant to meet at the 
same time, we suspect there is a deep reason 

/TeV^\' 

PA ~ -TIT- ■ (8-2) 



PI 



Indeed, p^'' is observationally about 2 meV, while TeV^/Mp; « 0.5 meV. 
Maybe that is why we see a coincidence M1061 . 
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Actually, an exact coincidence does not leave a window for structure forma- 
tion, which requires matter-dominated period. Fortunately, WIMP abundance is 
enhanced by weakness of the annihilation cross section, which goes like l/a^. 
This enhancement of matter relative to the triple coincidence gives us a window 
for matter domination and structure formation. May be that is why we seems to 
be in this triangle. But then, why is the baryon component also just a factor of 
five smaller than dark matter? Are they somehow related? 

Oh well, we know so little. 



9. Conclusions 

In my lectures, I tried to emphasize that we are approaching an exciting time to 
cross new threshold of rich physics at the TeV energy scale in the next few years 
at the LHC. At the same time, the dark matter of the universe is now established to 
be not made of particles we know, requiring physics beyond the standard model. 
The main paradigm for dark matter now is WIMPs, TeV-scale particles produced 
by the Big Bang which naturally give the correct order of magnitude for its abun- 
dance. Even though nature may be tricking us by this coincidence, many of us 
(including I) think that there is indeed a new particle (or many of them) waiting 
to be discovered at the LHC (or ILC) that tells us something about the dark side 
of the universe. If this is so, I would feel lucky to be born to this age. 



Appendix A. Gravitational Lensing 

Gravitational lensing is an important tool in many studies in cosmology and as- 
trophysics. In this appendix I introduce the deflection of light in a spherically 
symmetric gravitational field (Schwarzschild metric) 

Appendix A.l. Deflection Angle 
Using the Schwarzschild metric (c = 1) 

ds^ = dt^ — dr"^ - de^ - sin^ 6 dcf)"^ (A. 1) 

r r — rs 

where rs — 2Gjv™ is the Schwarzschild radius. The Hamilton-Jacobi equa- 
tioij^ for light in this metric is 



^"For an introduction to HamUton-Jacobi equations, 

|http: //hitoshi . berkeley ■ edu/22 lA/classical2 . pdf| 
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dSV r-rs (dSV 1 { dSV 1 { dSV 



r — rs \ dt J r \dr J r-^ \d9 J sin^ 6 \d(j> J 

^ 0. 

We separate the variables as 

S{t,r,0,cf,) = Slit) + S2{r) + SsiO) + Si{(l>) (A. 3) 

where 



r — rs \ dt J r \ dr J \ d9 J sin^ 9 \ del) J 

= 0. 

Because the equation does not contain t or (f) expHcitly, their functions must be 
constants, 

We can solve them immediately as 

Slit) = -Et, (A. 6) 

^4(0) - L,cj). (A. 7) 

Then Eq. dA. 2i becomes 

r — rs r \ dr ) r'^ \ m J sm 

The 6* dependence is only in the last two terms and hence 

is a constant which can be integrated explicitly if needed. Without a loss of 
generality, we can choose the coordinate system such that the orbit is on the x-y 
plane, and hence Lz — 0. In this case, Siicj)) — and 6*3(6') = L6. Finally, the 
equation reduces to 



r-rs fdS2, A' 



r — rs r \ dr 



^ - — ^ -jfM --T=0. (A. 9) 
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Therefore, 



= I \ r^^E^ - -r—^ dr. (A. 10) 

J \l (r- rsY r[r - rs) 

Since S{t,r,9,(t)) = 52 (r) ~ Et + L6, S2 can be regarded as Legendre transform 
52 (r, E, L) of the action, and hence the inverse Legendre transform gives 

2%^^^-<-- (A.n, 

oL 

Using the expression Eq. dA. lOK we find 

r Ldr 

0(r) = / , . (A. 12) 

Jr^ ^jE'^r^ - L^ir ~ rg) 

The closest approach is where the argument of the square root vanishes, 

E\t - L\,{r, - rs) - 0. (A. 13) 

It is useful to verify that the m = {rs = 0) limit makes sense. The clos- 
est approach is E^r"^ — L'^r^ = and hence r^ — L/E, which is the impact 
parameter. The orbit Eq. ( lA. 12\ is 



Ldr r rrdr 



Change the variable to r = rc cosh r], and we find 



(A. 14) 



e{r) = / '^^ '—!—— = / = 2arctantanh-^. (A. 15) 

Jg Tc cosh rjrc smn rj Jq cosh rj 2 

Hence tan | — tanh -I, and 



1 - tan^ 6'/2 _ 1 - tanh"^ ri/2 



rr 



cos 9 = = , = ; — = — . (A. 16) 

l+tan^6'/2 1 + tanh'' 77/2 cosh?] r 

Therefore rc — r cos 9 which is nothing but a straight line. 

To find the deflection angle, we only need to calculate the asymptotic angle 
9{r = 00). Going back to Eq. jA. 12b . we need to calculate 

9{oo) = / . (A. 17) 

Jr^ ^E^r^ - L^r{r - rs) 
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We would like to expand it up to the linear order in rs ^ re- If you naively 
expand the integrand in rs, the argument of the square root in the resulting ex- 
pression can be negative for r = rc < L/E. To avoid this problem, we change 
the variable to r = Vc/x: 

0(oo) = ^^^^"^ . (A. 18) 

Jo V E^r"^ — L'^rc{rc — rsx)x'^ 

Using Eq. ( lA. 13b . we write E'^r^ and obtain 

e{oo) = . (A. 19) 

Jo Vr2(l - .t2) - rcr5(l - x3) 

Expanding it to the linear order in rs / Vc, we find 

^ ' Jo \VT^ 2{l + x)VT^r, ^ ' ) 



2 r/ 



(A. 20) 



The deflection angle is A^? tt — 26{co) — 2^ = AGNtn/rc- It is easy to 
recover c = 1 by looking at the dimensions, and we find /S.9 ~ iGNm/c^Vc- 
It is also useful to know the closest approach Vc to the first order in m. We 

L 

E 



expand rc as rc = 4 + A. Then Eq. dA. 13l l gives 



A^^-2^^ + ^rs + 0{rsf ^Q, (A. 21) 

and hence 

rc = I - y + 0(rs)^ (A. 22) 

Appendix A.2. Amplification in Microlensing 

Once the deflection angle is known, it is easy to work out the amplification using 
simple geometric optics. Throughout the discussion, we keep only the first order 
in very small angles. Just by looking at the geometry in Fig. |20l the deflection 
angle is 

r — rn r — rn AG Mm 
AO = 01+62^ —-^ + —-^ ^ — —. (A. 23) 

di d2 r 
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Here, ro is the impact parameter. When rg = (exactly along the line of sight), 
the solution is simple: 



r(ro = 0) = i?o = W4GA,m 



did2 



(A. 24) 



This is what is called the Einstein radius, Rq in Paczynski's notation ll56l . For 
general rg, Eq. dA. 23b can be rewritten as 



r{r - ro) - ^ 0, 
which is Eq. (1) in the Packzynski's paper. It has two solutions 

r±{ro) = Uro±,fi+4R^o 



(A. 25) 



(A. 26) 



The solution with the positive sign is what is depicted in Fig.|20l while the solu- 
tion with the negative sign makes the light ray go below the lens. 



e, 



^____S- — — — > ' 






r 


T ^2 ' 

I 



)Star 



lens 



Fig. 20. The deflection of liglit due to a massive body close to tlie line of sight towards a star. 



To figure out the amplification due to the gravitational lensing, we consider 
the finite aperture of the telescope {i.e., the size of the mirror). We assume an 
infinitesimal circular aperture. From the point of view of the star, the finite aper- 
ture is an image on the deflection plane of size 6, namely the plane perpendicular 
to the straight line from the star to the telescope where the lens is. The ver- 
tical aperture changes the impact parameter rg to a range ro ± (5 (size of the 
mirror is (5 x [di + d2)/d2). Correspondingly, the image of the telescope is at 
r±{ro ±6) = r±{ro) ± Using the solution Eq. dA. 261 ), we find that the 

vertical aperture always appears squashed (see Fig. 12111. 



6 X 



dr 






= (5 X 


dro 






y/rg + 4i?g ± ro 
2v/^§+4lg 



< S. 



^'Note that this Taylor expansion is valid only when 5 ^ ro. For 5 ~ ro, we have to work it out 
more precisely; see next section. 
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On the other hand, the horizontal aperture is scaled as 

Sx—. (A. 27) 

Because the amount of light that goes into the mirror is proportional to the ellip- 
tical aperture from the point of view of the star that emits light isotropic ally, the 
magnification is given by 



A± = — 

To 



dr 



dro 



(Vrg+ 4j?g±ro)^ ^ 2rg+4fl§±2roVrg+4i?§ 



The total magnificiation sums two images, 

A^A,+A^^ ^3±3^ . 4±L (A. 28) 

with u = ro /-RoEl Basically, there is a significant amplification of the brightness 
of the star when the lens passes through the line of sight within the Einstein 
radius. 

Appendix A.3. MACHO search 

We estimate the frequency and duration of gravitational microlensing due to MA- 
CHOs in the galactic halo. The Large Magellanic Cloud is about 50kpc away 
from us, while we are about 8.5kpc away from the galactic center The flat ro- 
tation curve for the Milky Way galaxy is about 220 km/sec (see Fig. 6 in P8l ). 
The Einstein radius for a MACHO is calculated from Eq. ( lA. 24l l. 

= .l ^GNjn d,d, ^ ^ 2410^2 m ( ^\ ( ^) . (A. 29) 



To support the rotation speed of foo = 220 km/sec in the isothermal model 
of halo, we need the velocity dispersion cr = Voc /V2. The average velocity 
transverse to the line of sight is 

{vl+vl)=2<7'=vl. (A. 30) 

The time it takes a MACHO to traverse the Einstein radius is 

-^0 c c 1 n6 ( ""^ \ ^^^ ( \^did2\ .. , , , 

— = 5.6 X 10° sec — — -— , (A. 31) 

\MqJ V25kpcy' 



^•^The singular behavior for ro ^ is due to tlie invalid Taylor expansion in 5. This is practically 
not a concern because it is highly unlikely that a MACHO passes through with ro ^ <5o ■ Note that 
the true image is actually not quite elliptic but distorted in this case. 
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;±6; 



iro±6 




±6x/1>oV'o 



Fig. 21. The way the mirror of the telescope appears on the deflection plane from the point of view 
of the star. For the puipose of illustration, we took Rq = 2, ro = 3. 

about two months for m = Mq and di — d2 — 25 kpc. A microlensing event of 
duration shorter than a year can be in principle be seen@ 

The remaining question is the frequency of such microlensing events. It is the 
probability of a randomly moving MACHO coming within the Einstein radius 
of a star in the LMC^We will make a crude estimate. The flat rotation curves 

requires = ^ and hence the halo density p{r) — 4^^^° • '^^^ number 

2 

density of MACHOs, assuming they dominate the halo, is then n{r) = 4^g^^,.2 ■ 
Instead of dealing with the Boltzmann (Gaussian) distribution in velocities, we 
simplify the problem by assuming that v\ — v'^ + Vy — . From the transverse 
distance r\_ — ^ ^ y'^, only the fraction Ro/r± heads the right direction 
for the distance a At. Therefore the fraction of MACHOs that pass through the 
Einstein radius is 



We then integrate it over the depth with the number density. The distance from 
the solar system to the LMA is not the same as the distance from the galactic 

MACHO collaboration did even more patient scanning to look for microlensing events longer 
than a year ||67| . 




(A. 32) 
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center because of the relative angle a = 82°. The solar system is away from the 
galactic center by — 8.5 kpc. Along the line of sight to the LMA with depth 
R, the distance from the galactic center is given by ^ R^ + rl~ 2RrQ cos a 
with a = 82°. Therefore the halo density along the line of sight is 

'^^'^^ " 47rGjvm(i?2 + r| - 2RrQ cos a) 

The number of MACHOs passing through the line of sight towards a star in the 
LMA within the Einstein radius is 

dRn{r)2TTRoa 

Rlmc y2 

dR—- -— ^ -2TTRoaAt (A. 34) 

47rG'Arm(i?2 + r| - 2i?r0 cos a) 

Pakzynski evaluates the optical depth, but I'd rather estimate a quantity that is 
directly relevant to the experiment, namely the rate of the microlensing events. 
Just by taking At away, 

/•Rlmc 2 

rate = / dR — — f -2TrRoa 

Jq ATrGM'myR'^ + ~ 2i?r0 cos a) 



""'"^ vl RjRLMc-R) <y 

i?2 + r| _ 2RrQ cosay GnitiRi^mc c 

^R{Rlmc ~ R) dR 



c^GNmRLMC Jo R^ + r'^- 2RrQ cos a ' 



(A. 35) 



The integral can be evaluated numerically. For Rlmc — 50 kpc, tq ~ 8.5 kpc, 
a = 82°, Mathematica gives 3.05. Then with a = Woo/V2, ^'oo = 220 km/sec, 
we find 



rate = 1.69 x 10"'^ sec"' — = 5.34 x lO^Var"' — 

\ m J \ m 

Therefore, if we can monitor about a million stars, we may see 5 microlensing 
events for a solar mass MACHO per year, even more for ligher ones. 



Appendix A.4. Strong Lensing 

Even though it is not a part of this lecture, it is fun to see what happens when 
ro < 5. This can be studied easily with a slightly tilted coordinates in Fig.l22l 




Fig. 22. A slightly different coordinate system to work out the distortion of images. 



Using this coordinate system, we can draw a circle on the plane (x, y) = 
(xq, uo) + p{cos (j), sin (f)), and the corresponding image on the deflector plane is 
(i, y) = (io, j/o) + /5(cos 0, sin <j)) — j^^fj^ {x, y). The impact parameter is then 

To = ^JW^+ip which allows us to calculate r± (rp) using Eq. jA. 26\ for each (j). 
Obviously (j) is the same for the undistorted and distorted images. Fig. |23] shows 
a spectacular example with (a;o,?/o) = (IjO), -j^^^ ~ 1' P ~ ^-^^ Because 
p ^ ro, the Taylor expansion does not work, and the image is far from an ellipse. 




Fig. 23. A highly distorted image due to the gravitational lensing. Yellow circle is the undistorted 
image, while the two blue regions are the images distorted by the gravitational lensing. 



This kind of situation is not expected to occur for something as small as the 
mirror of a telescope, but may for something as big as a galaxy. When an image 
of a galaxy is distorted by a concentration of mass in the foreground, such as 
a cluster of galaxies, people have seen spectacular "strong lensing" effects, as 
shown in Fig.l24l 
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Fig. 24. A Hubble Space Telescope image of a gravitational lens fonned by the waiping of images of 
objects behind a massive concentration of dark matter. Warped images of the same blue background 
galaxy are seen in multiple places. The detailed analysis of lensing effects allows one to map out the 
mass distribution in the cluster that shows a smooth dark matter contribution not seen in the optical 
image. Taken from 1 107|. 
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